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LOWER LIMIT ON THE STRENGTH AND FILLING FACTOR OF EXTRAGALACTIC MAGNETIC FIELDS 
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ABSTRACT 

High energy photons from blazars can initiate electromagnetic pair cascades interacting with the 
cxtragalactic photon background. The charged component of such cascades is deflected and delayed 
by extragalactic magnetic fields (EGMF), reducing thereby the observed point-like flux and leading 
potentially to multi degree images in the GeV energy range. We calculate the fluence of 1ES 0229+200 
as seen by Fermi-LAT for different EGMF profiles using a Monte Carlo simulation for the cascade 
development. The non-observation of 1ES 0229+200 by Fermi-LAT suggests that the EGMF fills 
at least 60% of space with fields stronger than 0(1O -16 — 10~ 15 )G for life times of TeV activity 
of O(10 2 — 10 4 )yr. Thus the (non-) observation of GeV extensions around TeV blazars probes the 
EGMF in voids and puts strong constraints on the origin of EGMFs: Either EGMFs were generated 
in a space filling manner (e.g. primordially) or EGMFs produced locally (e.g. by galaxies) have to be 
efficiently transported to fill a significant volume fraction, as e.g. by galactic outflows. 
Subject headings: magnetic fields: origin - magnetic fields: cxtragalactic - gamma rays: galaxies - 
galaxies: active 



1. INTRODUCTION 

While magnetic fields are known to play a promi- 
nent role for the dynamics and in the energy budget 
of astrophysical systems on galactic and smaller scales , 
their role on larger scale s is still elusive (jWidrow 1 120021 : 
iKulsrud fc Zwc ibcl 200§). Extragalactic magnetic fields 
(EGMF) are notoriously difficult to measure and the 
data are incomplete. So far only in a few galaxy clus- 
ters observational constraints have been obtained, cither 
by observing their synchrotron radiation halos or by per- 
forming Faraday rotation measurements (RMs). Within 
galaxy clusters the inferred magnetic fields are between 
0.1 and 1.0 [iG on scales as large as 1 Mpc, and can be as 
strong as 30 /iG localized inside cluster cool cores. How- 
ever, as both observational methods need a prerequisite 
to measure magnetic fields (high thermal density for RMs 
and presence of relativistic particles for radio emission), 
they have been successfully applied only to high den- 
sity regions of collapsed objects as galaxies and galaxy 
clusters. Fields significantly below /iG level are barely 
detectable with these methods. Also other constraints, 
for instance the absence of distortions in the spectrum 
and the polarization properties of the cosmic microwave 
background radiation implies only a fairly large, global 
upper limit on the EGMF at the level of 10~ 9 G. 

The observed magnetic fields in galaxies and galaxy 
clusters are assumed to result from the amplification of 
much weaker seed fields. Such seeds could be created 
in the early universe, e.g. during phase transitions, and 
then amplified by plasma processes. Alternatively, an 
early population of starburst galaxies or AGN could have 
generated the seeds of the EGMFs at rcdshift between 
five and six, before galaxy clusters formed as gravita- 
tionally bound systems. In both cases, a large fraction 
of the material collapsing to today's visible structures 
will be seeded by such fields, but whether a significant 



fraction of the volume of the universe can be filled by 
such EGMFs depends on the efficiency of the transport 
processes involved. A quite different possibility is that 
the ejecta of AGN magnetized the intracluster medium 
only at low redshifts, and that thus the EGMF are con- 
fined within galaxy clusters and groups. The detection 
of EGMFs outside clusters is therefore crucial in discrim- 
inating different models for the origin of their seed fields. 

An alternative approach to obtain information about 
the EGMFs is to use its effect on the radiation from TeV 
gamma-ray sources. The multi- TeV 7-ray flux from dis- 
tant blazars is strongly attenuated by pair production 
on the infrared/optical extragalactic background light 
(EBL), initiating electromagnetic cascades in the inter- 
galactic space. The charged component of these cas- 
cades is deflected by the EGMF. Potentially observable 
effects of such electromagnetic cascades in the EGMF 
include the delayed "e choes" of multi- TeV 7-ray flares 
or gamma-ray bursts (|Plagal Il995t iMurase et all 12001 
and the appearance of extende d emission around ini- 
tially point-like 7-ray sources ([Aharonian et al.l 119941: 



Neronov fc Semikodl2007t iDolag et al.ll2009t lElviv et al.l 
20091: iNeronov et al.N2010al) . 



An additional way to derive lo wer limits on the EGMF 
has been pointed out recentl y bv INeronov fc VovkJ ff2?)10h 
and iTavecchio et all (|2010al ): Since the deflection of the 
cascade flux into an extended halo weakens the point- 
like image, the non-observation of TeV blazars in the 
GeV range by Fermi-LAT can been used to derive a 
lower limit on the EGMF. Particular suitable candidates 
are blazars with a very hard TeV spectrum like 1ES 
0229+200 that show a low intrin sic G eV emission. In 
this wa y, INeronov fc Vovkl (|2010f ) and ITavecchio et al.l 
(|2010aj) derived the lo wer bound B>5x 10~ 15 G on 
the EGMF. Note that lAndo fc Kusenkol (f20Toh found 
evidence for gamma-ray halos in the stacked images of 
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the 170 brightest active galactic nuclei observed with 
Fermi-LAT. The size of these halos are consistent with 
B ~ 10~ 15 G, however it is unclear if the excess found is 
real or due to the imperfect knowledge of the Fermi-LAT 
point sp read function (PSF) as argued bv lNeronov et all 
(|2010b[ ). All these analyses assumed a stationary source 
and did not discuss the effect of time delays induced by 
EGMFs. 

In this work, we will concentrate on 1ES 0229+200 and 
follow essentially the same assumptions about the source 
and th e sensitivity of Fermi-LAT as iTavecchio et al.l 
(|2010af ). We will improve on previous analyses in two 
respects: First, we use a Monte Carlo simulation for the 
development of electromagnetic cascades in the EBL that 
includes the effects of magnetic fields like synchrotron 
radiation and deflections of electrons. Moreover, we cal- 
culate the time delay that is induced by the EGMF de- 
flections. Second, we examine the influence of a more 
realistic, structured magnetic field on the EGMF limit. 

2. SIMULATION PROCEDURE 

We inject at the rcdshift z — 0.14 photons with an 
energy distribution given by J- oc E~ 2 ^ 3 and a maximal 
energy i? max = 20TeV. Such a hard injection spectrum 
gives a reasonable fit to the HESS data an d is consistent 
with recent SWIFT dBurrows etafl 120041 ) observations 
as arg ued by ITavecchio et al.1 (|2009l ). We assume also a 



relatively low Lorentz factor, T = 10, i.e. 0j C t = 6 
order to derive a conservative lower bound on the EGMF. 

We follow electromagnetic cascades initiated by these 
high ene rgy photons with the M onte Carlo code intro- 
duced by |KachehieIiEaD j2009|). Since Fermi-LAT pro- 
vides only an upper limit on the photon flux from 1ES 
0229+200 it is sufficient to assume a radial-symmetric 
halo and to calc ulate the deflectio ns and time delays 
as described in (|Dolag et al.l 120091) . We use as EBL 
background the best-fi t model from the calculations of 
IKneiske fc Doll (pfTTol) . 

We divide the photons arriving at the observer into 
a point-like flux and a halo outside a circle with angle 
t9* around the source. For Fermi-LAT, we use for i9* 
an analytical approximation of the angle 1^95 containing 
95% of the flux emitted by a point-source given by $95 « 
1.68°(i;/GcV)- - 77 + 0.2°exp(-10GeV/£:). Above 300 
GeV, t?95 is set equal to 0.11°, a value typical for the 
resolution of atmospheric Chercnkov telescopes. 

3. RESULTS FOR DIFFERENT EGMF MODELS 

In Figs. [T] and [2] we show our results for the fluence con- 
tained inside the 95% confidence contour of the PSF of 
Fermi-LAT. Additionally, these figur es contain the HESS 
observations (jAharonian et al.l 120071 ) as black dots with 
error bars and t he Ferm i-LAT upper limits derived by 
ITavecchio et al.l ()2010al ). The fluences have been nor- 
malized fitting them to the HESS data. 

In Fig. [TJ we have used a uniform magnetic field so 
that our results can be directly compared w ith the an- 
alytical estimates of ITavecchio et all (|2010al ). Note that 
a turbulent field with correlation length L cl - much larger 
than the mean free path lie of electrons in the Thomson 
regime, L CI lie ~ 1 kpc, is well approximated by a 
uniform field. For smaller correlation lengths, L CY <C he, 
the electron diffuses in the small-angle deflection regime, 
requiring larger magnetic fields for the same deflection 
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Fig. 1. — Fluoncc contained inside the 95% confidence contour of 
the PSF of Fermi-LAT as function of energy together with Fermi- 
LAT upper limits and HESS observations for a uniform EGMF with 
strengths varying from B = 10" 16 G to B = 10" 13 G with £ max = 
20TeV (solid) and 100 TeV (dashed). The direct component for 
B = 10 -14 G is also shown. 
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Fig. 2. — Fluence contained inside the 95% confidence contour 
of the PSF of Fermi-LAT as function of energy for a EGMF with 
top-hat profile and filling factor / varying from / = 0.1 to / = 0.9 
with _B max = 20 TeV (solid) and 100 TeV (dashed). 

angle. Demanding that the cascade flux is below the up- 
per limits of Fermi-LAT leads to a lower limit on the 
magnetic field strength of <~ 10~ 14 G. For this case the 
direct component, i.e. photons arriving at the detector 
without cascading, is also shown. Note that for small 
-E-max the transition from the direct to the cascade con- 
tribution leads to a break at ~ TeV in the spectrum, as 
suggested by the HESS data. 

While our li mit agrees reasonably we ll with the analyt- 
ical estimate of ITavecchio et all (|2009D . the shape of the 
cascade flux obtained differs. There ar e several reasons 
respon sible for these differences: First, ITavecchio et al.l 
(|2009l) assume that the spectral shape of the cascade 
flux below the threshold energy ~ 10 11 eV is given 
by J 7 oc E~ - 5 for negligible magnetic fields. Such a 
slope typical for the regime of Thomson cooling is how- 
ever restricted to energies E < 10 8 cV, while at higher 
energies a plateau T oc E a with a ~ —0.9 is ex- 
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pected ()Berezinskv et al1l2010l) . Second, deflections in 
the EGMF lead even for isotropic emission to extended 
images o f point-like sources. T his effect has been ne- 
glected in l^ecchioeEaD J2QQ^). Third, using full prob- 
ability distributions for the interactions there is a non- 
negligible probability for photons not to interact, espe- 
cially towards the low-energy end of the injected energy 
range. Finally, the energy dependent PSF of Fermi in- 
troduces an artificial energy dependence of the point-like 
flux T cx . Note that an increase in -E max from 20 to 
100 TeV reduces the limit on the magnetic field strength 
to ~ 5 x 10" 15 G, see Fig. 1, while a further increase of 
fi mas strengthens the limit again. The counter-intuitive 
behavior between 20 to 100 TeV is caused by the dom- 
inance of direct photons at the high-energy part of the 
spectrum for small E max . 

We discuss next the consequence of time delays in- 
duced by the EGMF. We stress that the non-observation 
of 1ES 0229+200 by Fermi-LAT is a signature for non- 
zero EGMFs, even if the delay of the Fermi signal would 
exceed the life time r of 1ES 0229+200 as TeV source. 
However, the numerical value of the limit deduced using 
deflections would be modified. For B = 5 x 10" 15 G, we 
find as average time delay 1 x 10 6 , 1 x 10 5 , and 5000 yr for 
the three energy bins (0.1-10), (1-10), and (10-30) GeV. 
Thus the high-energy bin of the Fermi data can be used 
to derive the constraint as above, if 1ES 0229+200 is 
a relatively stable source over a time scale longer than 
10 4 yr. Even if r is much smaller, say around 100 yr, 
the lower limit on B weakens only by a factor 10 to 
B = £>(10" 16 G). 

Since the EGMF is strongly structured, one may won- 
der how a non- uniform field modifies this limit. In par- 
ticular, we want to address the question whether the 
presence of relatively strong fields concentrated inside 
cosmic structures like filaments could mimic the effect 
of an EGMF present also in voids. As simplest possi- 
ble test, we use first a top-hat profile for the structure of 
the EGMF: We set the field strength to zero in a fraction 

I — / of space and use a value which in general is assumed 
to be representative for filaments, B = 10~ 10 G, in the 
remaining part. For the separation of the peaks we use 
D = lOMpc motivated by the typical distances between 
cosmological structures, although the exact value of D 
plays no role as long as (1 — f)D <C / 7 , with Z 7 as mean 
free path of photons. The dependence of the fluence con- 
tained inside the PSF of Fermi-LAT on the filling factor 
/ is shown in Fig. [2] To be consistent with the Fermi up- 
per limits, sufficiently strong magnetic fields should fill 

> 80% of space. The derived limit on the filling factor 
is practically independent of the source life time r and 
B, as long as the field is stronger than >5xlO" 15 G. As 
in the previous case, by assuming a higher injected E max 
the required filling factor is slightly reduced to 60% . 

The failure of strong fields filling only a small frac- 
tion of the universe to suppress sufficiently the point-like 
cascade flux can be understood as follows: The HESS ob- 
servations of 1ES 0229+200 cover the energy range 0.5- 

II TeV. In the same energy range, the mean free path 
Z 7 of VHE 7-rays through the EBL varies between 1000 
and 50 Mpc and is thus always much larger than the typ- 
ical extension of regions with large fields, (1 — f)D. For 
the energies considered, the cascade consists typically of 
only three steps, 7 — > e ± — > 7. Since the mean free path 
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Fig. 3. — The EGMF component perpendicular to the line-of- 
sight towards 1ES 0229+200 from four different MHD simulations. 
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Fig. 4. — Cumulative volume filling factor C(B) for the four 
different EGMF models found in MHD simulations. 



he of electrons in the Thomson regime is very small, 
lie ~ 1 kpc, all cascades with electrons created outside 
the strong-field regions are undeflected. Thus it is not 
possible to trade smaller values of / against larger val- 
ues of B: Increasing the field strength beyond ~ 10" 13 G 
leads only to an increase of the deflection, while the frac- 
tion of cascades deflected outside the Fermi PSF remains 
constant. 

To mimic the even more complex situation in pres- 
ence of cosmological structures, we have used magnetic 
fields derived from several cosmological MHD simulations 
using constrained initial conditions (Dolag et al. 2004, 
Donnert et al. 2009; for alternative simulations see e.g. 
Ryu et al. 2008). Therefore they reproduce fairly well 
the true large scale matter distribution, out to a distance 
of w 114 Mpc, cf. for more details Dolag et al. (2005). 
We extracted the EGMF as predicted by these simula- 
tions along the line of sight towards the position of 1ES 
0229+200 for several models for the magnetic seed field. 
A comparison of the EGMF component perpendicular to 
the line-of-sight towards 1ES 0229+200 obtained in the 
models is shown in Fig. [3] 

Most magnetic seed models used in cosmological simu- 
lations assumed a space filling (e.g. primordial) magnetic 
seed field of order 10 _12 G, which inside galaxy clusters is 
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Fig. 5. — Fluence contained inside the 95% confidence contour 
of the PSF of Fermi-LAT as function of energy for EGMFs from 
four different MHD simulations with -E m ax = 20TeV (solid) and 
100 TeV (dashed). 

amplified by the collapse and additionally by shear tur- 
bulence up to the observed fiG level. For the purpose 
of this work we have selected the available simulation 
with the lowest primordial seed field of 10 -13 G (MHDx), 
which inside voids is even diluted to values of order of 
10~ 14 G due to cosmic expansion. To construct a sec- 
ond model with magnetic fields as low as 10 -17 G inside 
voids, we scaled down the magnetic field of the original 
simulation by a factor proportional to the local density 
to the power 2/3 (Model 3x) with the fix point at the 
core density of galaxy clusters. Thereby, the Faraday 
RMs of the newly constructed model within galaxy clus- 
ters does not significantly change, whereas the magnetic 
field in filaments and voids drop by nearly three orders 
of magnitude compared to the original simulation. 

Models which intrinsically produce lower magnetic 
fields (and filling factors) in filaments and voids are those 
where the seed field originates from the galaxies formed. 
In such simulations, galaxies are identified at early times, 
and their halos are magnetized according to the magnetic 
field observed in nearby galaxies. So far such simulations 
do not include any transport processes for the magnetic 
field (like galactic outflows or AGN activity) beyond the 
large scale potential well in which galaxies live. Never- 
theless, these models reproduce well the observed mag- 
netic field properties within galaxy clusters, since the gas 
component of the galactic halos is stripped during the 
formation process of the cluster and the magnetic field 
gets well mixed and subsequently amplified by the turbu- 
lence within the intercluster medium. The result of such 
a simulation is the model (MHD Gal3) from Donnert et 
al. 2009. As galaxies in filaments form later than their 
counterparts in high density environment, we have also 
extracted the magnetic field from a model which subse- 
quently put magnetic fields into newly formed halos also 
at lower redshift (MHD Gal5). The main difference to 
the previous model is the presence of more magnetized 
structures within filaments whereas the magnetic field 
properties in galaxy clusters do not change significantly. 

The simulation used covers only up to distances of 
114 Mpc, and therefore we are forced to extend the mag- 
netic field profile beyond that point. For this purpose we 



assume that the extended field has the same statistical 
properties as the simulated patch. In practice, we simply 
mirror and repeat the field beyond the first 114 Mpc. In 
Fig. [4] we show the cumulative volume filling factor, i.e. 
the fraction C(B) of volume occupied by regions with 
field strength > B. The vertical and horizontal bands 
stand for the minimum magnetic field and filling factor, 
respectively, required to avoid the Fermi limit, as ob- 
tained in Figs.[T]and[2j 

The fluence contained inside the PSF of Fermi-LAT for 
the EGMF models MHDx, MHD Gal3 and MHD Gal5 is 
shown in Fig. [5] From the results for the top-hat model, 
we anticipate that the EGMF in the models MHD3x and 
MHD Gal3 is too weak and does not lead to a signif- 
icant modification of the fluence, see Fig. [4j Since the 
model MHD Gal5 is not far from the intersection of the 
two bands, we expect it not to lie far from the Fermi 
limit. In Fig. [5] it is seen how in the particular case of 
i? mM = 100 TeV this model is marginally consistent with 
the Fermi-LAT observations. 

Finally, the fluence is several orders of magnitude be- 
low the observational limits for fields as strong as in the 
model MHDx. In such a case, extended image s in the 
GeV range as predicted e.g. by INeronov et al.l (|2010aD 
will be impossible to observe with Fermi-LAT. 

We would like to point out, that on the other hand, all 
these models are in agreement with observations of mag- 
netic fields in galaxy clusters - which so far is the only 
place where the strength and structure of magnetic fields 
has been directly detected - and therefore demonstrate 
the constraining power of using high energy photons to 
infer the origin and the processes leading tho EGMFs. 

4. SUMMARY 

We have calculated the fluence of 1ES 0229+200 as 
seen by Fermi-LAT using a Monte Carlo simulation for 
the cascade development. We have discussed the effect of 
different EGMF profiles on the resulting suppression of 
the point-like flux seen by Fermi-LAT. Since the electron 
cooling length is much smaller than the mean free path of 
the TeV photons, a sufficient suppression of the point-like 
flux requires that the EGMF fills a large fraction along 
the line-of-sight towards 1ES 0229+200, / > 0.6. The 
lower limit on the magnetic field strength in this volume 
is B ~ 0(1O- 15 )G, assuming 1ES 0229+200 is stable 
at least for 10 yr, weakening by a factor of 10 for r = 
10 2 yr. These limits put very stringent constraints on the 
origin of EGMFs: Either the seeds for EGMFs have to 
be produced by a volume filling process (e.g. primordial) 
or very efficient transport processes have to be present 
which redistribute magnetic fields that were generated 
locally (e.g. in galaxies) into filaments and voids with a 
significant volume filling factor. 
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Note added. After the completion of this work the 
preprint of Tavecchio et al. (2010b) appeared in which 
these authors study the four BL Lac objects (RGB 
J0152+017, 1ES 0229+200, 1ES 0347-121 and PKS 0548- 
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322) and derive a bound of B ~ 5 x 10 15 G on a space filling magnetic field. 
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